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Abstract 

We study the hadronic spectrum in quenched lattice QCD using 
the improved Wilson fermion action proposed in [y at /3 = 5.7 and 
(3 = 6.0. We find a systematic reduction of the finite spacing effects 
compared to the results obtained by using the standard Wilson action. 



1 Introduction 

Most of the recent results about quenched QCD spectroscopy g] exhibit 
the scahng behaviour expected near the continuum hmit 0. However, devia- 
tions are still observed in the phenomenologically relevant heavy quark mass 
region, and in these cases progress using the standard Wilson action seems to 
be slow. Moreover, an important self-consistency check of lattice computa- 
tions is the agreement between results coming from different discretizations. 

For the above two reasons, we decided to study lattice spectroscopy by 
using an improved version of the Wilson fermion action |jl[, inspired by 
Symanzik's perturbative procedure of eliminating finite-a corrections @]. 
Other proposals along the same line are the ones of ref. |^ and ref. 0, 
which, by combining the improved action with a suitable modification of the 
fermion operators, eliminates all terms of order a in the hadronic matrix el- 
ements. Some results obtained by using the actions proposed in refs. [|5|, [^ 
are already available and we remand to the references 0], [§], @ for their 
discussion. 

There are two points to be considered, in relation to an improvement 
program based on perturbative arguments. The first one is the real source of 
the scaling violations, which can be completely non-perturbative, the other 
one is the practical applicability. Indeed, the success of an improved action 
depends on the interplay between the reduction of 0{a) effects, and the 
possibly induced new systematics, and it is not guaranteed a priori. To 
understand the real effectiveness of the action |l[ we study in this paper 
/3 = 5.7 and (3 = 6.0. From the results for the standard Wilson case we know 
that the finite spacing effects are strong at /5 = 5.7, while /3 = 6.0 seems to 
be the onset of scaling for many relevant quantities. We thus use the data at 
P = 5.7 to demonstrate the suppression of the most evident scaling violations 
which occur using the Wilson action at the same value of the coupling, and 
the ones at /5 = 6.0 to explore the possibility to actually improve the current 
Wilson results. 

We introduce the action we use, and summarize its basic properties in the 
next section. Then, we describe the numerical simulation and the analysis 
procedure. Finally we discuss the results for the hadron spectrum and the 
meson decay constants. In the following we keep as a reference the results 
obtained with the standard Wilson action at /5 = 5.7, /3 = 6.0,/? = 6.3 of 



T0| , |Tl|. Our conclusions are based on this comparative analysis. 



2 The action 

We use the improved ferinion action proposed by Hamber and Wu [^. A 
next nearest neighbor term is added to the Wilson fermion action: 

S = Sg + Sw + Sii (1) 

Sg and Sw are the usual pure gauge and Wilson fermionic term, while the 
new term reads: 

Sii = Yl ^ri{C - D-ff,)Un,i^Un+^,,^,iJn+2fi + h.c. (2) 

The choice C = —1/Akr cancels at tree level the term 0{p^) in the inverse 
propagator Sp^{p), D is free. We choose D = —l/8kr which cancels also the 
terms 0{p^) in the inverse fermion propagator, and set r = 1. In this way 
the fermionic part of the action coincides with the Eguchi-Kawamoto one. 

We refer to the original references |]l|, and to ref. for a detailed dis- 
cussion. It is however interesting pointing out here that this action exhibits 
positivity violations stronger that the ones of the standard Wilson case. The 
zero- momentum free quark propagator Sf{po, 0) has four poles for any k 
value, two of them complex, the other two turning from complex to real for 
k >~ 0.14.( The critical value kc , defined by mq{kc) = 0, ruq being the 
quark mass, is 1/6.) So, the free propagator, computed by Fourier trans- 
forming Sf{p) to the coordinate space, deviates from a simple exponential 
behaviour. It shows a clear ripple at small /c's (i.e., in the region of four com- 
plex poles), which is more evident in small lattices. The effect progressively 
disappears by increasing the size of the lattice, and approaching kc, suggest- 
ing that no problem arises in the continuum limit. Indeed, in the continuum 
limit mpHYsO' = it^latt -^ 0, while k -^ k^.. So, in order to define properly 
the a — i> limit only one real pole (the one corresponding to the physical 
mass) should be zero as k ^ k^, while the other real pole, and the real parts 
of the complex ones, should be different from zero. It is trivial to show that 
this is the case. Turning to the amplitudes, no problem should arise with 
the ones involving physical status: violations of positivity only affects the 
residuals of non-physical poles. 

The above informal discussion, which leads to the conclusion that the 
spurious poles are harmless, concerns the free case. However, the results for 
the standard Wilson action suggest that what has been found in the free case 



is general: for the Wilson action, the free propagator has a complex pole for 
K > 1/6, and a real one otherwise, and it has been shown |T^ that the same 



pattern of violation of positivity is maintained also in the interacting case. 
In conclusion, a more detailed investigation of the positivity properties of the 
improved actions would be welcome, but at the moment we have no reason 
to think that their continuum limit is suspect. 

Anyway — even if this may seem rather paradoxical — an action sup- 
posed to reduce finite size effects potentially leads to a peculiar finite spacing 
systematics, whose actual impact is not possible to estimate a priori. Hap- 
pily enough, positivity violations affected the results only at /3 = 5.7, for the 
L = 1 mesons. 

3 The Numerical Simulation and the Data 
Analysis 

We used the Wilson pure-gauge action, and we have thermalized the 
gauge fields at /5 = 5.7 and (3 = 6.0. We used a standard Metropohs-5 
hits algorithm to generate the background gauge field configurations. The 
onset of the scaling region for the standard Wilson action seems to be around 
P = 6.0 while the results at /5 = 5.7 are definitively affected by finite spacing 
effects. Thus, the values we choose are good to assess the effectiveness of the 
method. 

We explored the heavy quark region, which is appropriate to test this 
action. We choose four equispaced K, ranging from 0.186 to 0.198 at (3 = 5.7, 
and from 0.176 to 0.188 at j3 = 6.0. This gives the ratio vr^/p^ between 0.4 
and 0.8 at /5 = 5.7 and between 0.6 and 0.9 at /3 = 6.0. 

One configuration every 800 iterations was sampled for the propagator 
inversion. At (3 = 5.7 we have computed (140 x 4 quark masses) propagators 
on a 24 X 12^ lattice, at (3 = 6.0 we have (160 x 4) propagators on a 32 x 12^ 
lattice. We also collected (180 x 3) propagators on a 18 x 9^ lattice (at the 
three bigger masses) to monitor the finite size effects. 

All the measures were performed according to the methods introduced 
TU| and reviewed in []T3[, which we briefly summarize in the following. 



m 



We compute the fermionic propagators by a preconditioned minimal residue 
algorithm, as modified for parallel updating. We apply the incomplete LU 



decomposition as a preconditioner, making use in this step of the standard 
Wilson operator. In practice the original equation for the improved propa- 
gator 

DX = B (3) 

is replaced by 

U-^L-^DX = U-^L-^B (4) 



where L and U are defined in the usual way |T0[ in term of the Wilson 
fermionic operator. We have indeed found also with the improved action a 
considerable speed-up in the convergence rate by using this preconditioning. 

We use smearing |jl3| [|l^: after the necessary gauge fixing we solve the 
Dirac equation with an extended source (a 3^ spatial cube). The propagators 
computed in this way can be smeared on the final point too. From the quark 
propagators we form the contractions needed to evaluate the hadronic Green 
functions Gq (smeared only at the origin) and G^ (at the origin and at the 
final point). We preferred the second ones {Gx, smeared at the origin and 
in x) since Go sometimes tends to amplify the problems related to the lack 
of positivity of this action, but the difference is small and on the overall the 
results obtained from Gq and G^ are fully consistent. A joint fit of Gq and G^ 
does not further reduce the errors, which are always computed by a standard 
jacknife analysis. 

The results for the hadron masses quoted in Table || and in Table ^ are 
from a single mass fit (for t > 5) of the Green functions smeared both at the 
origin and at the final points. We systematically controlled the stability of 
the results by performing two particle fits as well (in this case we included 
points starting form t =~ 3), and we cross checked with the local fits. As 
an example of the overall quality of the data we show in Fig. 1 the effective 
masses of the proton at /3 = 5.7 and in Figs. 2 the fits at /5 = 6.0 for the 
proton and the pion. 

Using the data on the 9^ x 18 lattice we checked that finite size effects are 
small. The big symbols in Fig. 1 which perfectly circle the results obtained 
on the 12^ x 24 lattice are from the smaller one. The full set of results on 
the 9^ X 18 lattice is shown in Table ^j. (We note that the plateau in the 
effective masses in this lattice is not as extended as in the bigger one, so 
in this case we had to rely on two particle fits.) The results reported in 
Table ^ agree with the ones of Table |l|, thus demonstrating the good control 
we have over finite size effects. We are confident that, since the time extents 





k 


Particle 


0.186 


0.190 0.194 0.198 


TT 


1.0267( 28) 


0.8661( 29) 0.6975( 31) 0.5069( 40) 


7f 2 


1.0239( 39) 


0.8624( 42) 0.6925( 50) 0.4975( 73) 


P 


1.1432( 40) 


1.0109( 44) 0.8851( 59) 0.7676(132) 


P 


1.1451( 44) 


1.0141( 51) 0.8905( 73) 0.7811(161) 


P 


1.7972(108) 


1.5811(116) 1.3616(141) 1.1186(277) 


A 


1.8646(138) 


1.6665(150) 1.4761(198) 1.2990(557) 



Table 1: Hadron masses as a function of /c at /9 = 5.7 on the 24 x 12^ lattice 
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Particle 


0.186 0.190 0.194 


71 
TT 
P 
P 
P 

A 


1.0199( 30) 0.8638( 31) 0.6970( 37) 
1.0207( 40) 0.8623( 46) 0.6714( 59) 
1.1326( 43) 1.0073( 47) 0.8863( 57) 
1.1327( 48) 1.0052( 56) 0.8814( 75) 
1.6807(182) 1.5465(126) 1.3505(100) 
1.6541(294) 1.5798(219) 1.4707(212) 



Table 2: Hadron masses as a function oi k a.t (3 = 5.7 on the 18 x 9^ lattice 



in physical units of the lattices at /3 = 5.7 and /3 = 6.0 are roughly the same, 
the masses evaluated at /3 = 6.0 do not suffer from finite size effects as well. 
(However, the number of points in the space directions is the same for the 
two lattices, so some problem at least with the smallest mass at /5 = 6.0 
cannot be excluded.) 

At f3 = 5.7 we do not quote the results for the al and bl particles. 
These mesons deserve some separate comments, since their behaviour turned 
out to be quite peculiar. Their Green functions can even change sign, thus 
indicating violations of positivity ( amplitudes of different signs in different 
channels, or manifestations of the spurious complex poles). As a consequence, 
the effective masses have a wiggle which becomes less and less apparent as 
A; —> fee or a ^ 0, in agreement with the qualitative discussion of Sect. 2 above. 
In fact, the results at /5 = 6.0 for the fits of the L = 1 mesons are reasonable. 
(Of course in this case we do not have any hint about residual finite size 
effects, and it is clear that having amplitudes of comparable magnitude and 





k 


Particle 


0.176 


0.180 0.184 


0.188 


71 


1.0014( 15) 


0.8179( 17) 0.6287( 23) 


0.4204( 39) 


n 


0.9976( 17) 


0.8141( 21) 0.6252( 29) 


0.4170( 52) 


P 


1.0424( 19) 


0.8712( 24) 0.7024( 33) 


0.5333( 63) 


P 


1.0408( 20) 


0.8691( 25) 0.6990( 36) 


0.5238( 76) 


al 


1.2794(142) 


1.1203(187) 0.9676(285) 


0.8150(570) 


hi 


1.2772(140) 


1.1159(177) 0.9553(262) 


0.7873(557) 


P 


1.6356( 43) 


1.3734( 55) 1.1126( 80) 


0.8458(150) 


A 


1.6602( 48) 


1.4061( 61) 1.1584( 89) 


0.9120(183) 



Table 3: Hadron masses as a function of A; at /3 = 6.0 on the 32 x 12^ lattice 



opposite signs is the worst possible situation from this point of view.) As 
a general comment, the orbitally excited hadrons are difficult to treat, also 
with the ordinary Wilson action, and their signal is intrinsically noisy: a 
recent discussion, together with new results, can be found in ref. []T4| . 

We conclude the description of the analysis procedure by discussing the 
computations of the amplitudes of the Green functions in the fundamental 
channels, which are important for the estimates of the decay constants. We 
are interested in the local amplitudes, which, as discussed in ref [|ll| can 



be recovered from the ratio of Gq to Gx in the hypothesis of long-distance 
factorization of the Green functions. We checked that the amplitude of the 
fit of the ratio has a smaller statistical error compared to the ratio of the 
amplitudes of the fits, typically by a factor 2. This is understandable, because 
by fitting the ratio we get rid of the coherent fiuctuations. So the errors we 
quote are obtained in this way (the central values are basically the same). 
At P = 5.7 we also computed 24 propagators with a point source, so we 
can evaluate directly the local amplitudes which turn out to be in complete 
agreement with the ones coming from the fit of the ratios. The summary of 
the results for the amplitudes is given in Tables ^ and |[ 





A; 


Particle 


0.186 0.190 0.194 0.198 


TT^ 


2.140( 57) 1.950( 53) 1.834( 56) 1.876( 78) 


7r2 


2.123(100) 1.904(105) 1.739(112) 1.647(141) 


Tfl 


0.308( 14) 0.202( 10) 0.119( 7) 0.058( 5) 


7f2 


0.299( 26) 0.192( 21) 0.112( 16) 0.058( 14) 


P' 


0.977( 38) 0.779( 33) 0.611( 33) 0.484( 53) 


2 


0.912( 61) 0.714( 51) 0.538( 42) 0.358( 54) 



Table 4: Amplitudes in the fundamental channel as a function oik dX j3 
on the 24 x 12^ lattice, from the ratios' fit^ and from local operators^ 
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k 


Particle 


0.176 


0.180 0.184 


0.188 


TT 
TT 

P 


0.675( 15) 
0.183( 4) 
0.340( 8) 


0.540( 14) 0.427( 13) 
0.116( 3) 0.065( 2) 
0.242( 7) 0.160( 6) 


0.342( 17) 
0.027( 1) 
0.092( 6) 



Table 5: Amplitudes in the fundamental channel as a function oi k aX (3 = 6.0 

4 Hadron masses, decay constants and scal- 
ing behaviour 

We begin the discussion of our results with the chiral behaviour. Since 
we were particularly concerned with heavy masses, the chiral extrapolation, 
especially at /? = 6.0, is somewhat delicate. For heavy flavors, the data 
should follow the predictions of the potential models for quarkonium, whose 
curvatures as a function of the quark mass are different from the ones in 
the chiral limit. Consequently, it is difficult to decide when a fit, even if 
satisfactory on statistical grounds, correctly describes the chiral behaviour. 
With this warning in mind, we fitted the squared hadron masses as second 
order polynomials in the quark masses, defined as rriq = 2/3(1/ k — 1/kc). 
This turned out to be the most effective parametrization, in agreement with 
previous experiences with heavy masses. 

The first step is the determination of kc which is obtained by demanding 



Particle 


2 


Sl S2 


vr 





2.687( 49) 2.492( 84) 


TT 


0.000( 0) 


2.594( 78) 2.747( 219) 


P 


0.387( 38) 


1.984( 278) 3.356( 517) 


P 


0.422( 48) 


1.775( 358) 3.716( 663) 


P 


0.551( 114) 


7.648( 828) 3.635(1572) 


A 


1.150( 282) 


5.377(2224) 7.322(4247) 



Table 6: Coefficients of the extrapolating polynomial for the hadron masses 
at (3 = 5.7 



m^{kc) = 0. We get at /5 = 5.7 



m. 



2.687(49)mg + 2.492(84)m; 



and kc = 0.20333(12) At /3 = 6.0 we have 



m. 



2.032(40)mg + 3.504(46)m; 



(5) 



(6) 



and kc = 0.19216(12) 

It is interesting to compare with the analytical predictions for kc which 
are poorly reproduced on the lattice in the standard Wilson case. We have 
a 9 percent deviation at /5 = 6.0 and 15 percent at /? = 5.7 from the ana- 
lytical result kc = 1/6(1 + 1/18^^^ + 0{g^)) (second entry of ref. |jl|), to be 
compared with 12 percent and 17 percent for the standard Wilson action at 
the corresponding /3 values. 

The results for the pion at (3 = 5.7 and (3 = 6.0 with the results ffis 
superimposed are shown in Fig. 3. In the same figure we also show the 
results for the vf's, which turn out to be in full agreement with the ones for 
the IT. 

Once kc is determined, we go ahead with ffis in other channels. As al- 
ready said, they are completely satisfactory on statistical grounds, and do 
not deserve any special comments, so we simply quote in Table ^ and the 
coefficients of the extrapolating polynomials. 



hadron 



Mq + sinig + S2im?q 



(7) 



for the two (3's. 



Particle 


2 


Sl 


S2 


vr 





2.032( 40) 


3.504( 46) 


TT 


0.000( 2) 


1.996( 31) 


3.541( 85) 


P 


0.119( 10) 


1.872( 63) 


3.658( 108) 


P 


0.103( 12) 


1.971( 78) 


3.474( 136) 


al 


0.436( 146) 


2.742( 942) 


3.217(1698) 


61 


0.362( 136) 


3.180( 864) 


2.520(1548) 


P 


0.290( 36) 


4.926( 198) 


8.047( 323) 


A 


0.422( 51) 


4.718( 316) 


8.194( 541) 



Table 7: Coefficients of the extrapolating polynomial for the hadron masses 
at /3 = 6.0. 



We remark again the agreement between the results for the vr and the n 
particles, and for the p and the p particles. 

We note that at /5 = 5.7 the result mp(0) for the extrapolated p mass 
lies in between the Kogut-Susskind WB result and the Wilson one, support- 



ing the intuitive picture of an "effective spacing" between a and 2a for the 
improved action: 

mp{0){W) = 0.185(18) < mp{0){WI) = 0.387(38) < mp{0){KS) = 0.76(5) 

(8) 
At (3 = 6.0 the results for Wilson and Wilson improved are mutually consis- 
tent, while the result for Kogut-Susskind is getting closer to them (the ratio 
of the K-S resuhs to the Wilson one is 1.14(7) at /5 = 6.0): 

mp{0){W) = 0.111(5) ~ mp{0){WI) = 0.119(10) < mp{0){KS) = 0.144(27) 

(9) 
As already noticed in [16|, (3 = 6.0 seems the onset of the region in which 



the details of lattice discretization are forgotten. 

The comparison of the results obtained with the Wilson action and with 
the improved one is better done by using adimensional quantities: for this 
purpose we give in Table | the values of some relevant ratios. For the string 
tension at /3 = 5.7 we use the result of ref [0: aa^ = 0.056(2). At /? = 6.0 
and (3 = 6.3 we get, by analyzing the same configurations we were using in 
|n|, §, aa^ = 0.0471(35) and aa^ = 0.01704(8), respectively 0. (Let us 
remind that the experimental values to be compared with are p/p = 1.22, 
p/a = 2.24, p/a = 1.83, assuming the accepted estimate for a of 420 MeV.) 



10 





6.0 5.7 


6.3 W 6.0 W 5.7 W 


v/p 

p/a 


1.561(117) 1.193(136) 
2.481(179) 2.124(220) 
1.590( 89) 1.780( 88) 


1.247( 95) 1.298( 53) 1.451( 75) 
2.220(166) 1.993(101) 2.217( 93) 
1.780( 64) 1.535( 66) 1.528( 49) 



Table 8: Relationship between the extrapolated values for mg = of the p 
and proton masses and the string tension 

It is remarkable the full agreement of the improved results at /3 = 5.7 
with the ones obtained with the standard Wilson action at /5 = 6.3, which 
in turn match the experimental results. At /3 = 6.0 the situation is less clear 
(the p is the same as the normal Wilson case at the same (3, the proton is 
closer to the one at (3 = 6.3, but as previously discussed we do not have a 
safe estimate of the systematic errors induced by the chiral extrapolation at 
/3 = 6.0) 

In the following we will discuss finite mass data, using occasionally the 
extrapolated values for the p mass to set a common scale for the results 
obtained at different /3 values. 

An interesting quantity to look at is the splitting between pseudoscalar 
and vector mesons with the same flavour content. Experimentally, My — Mpg 
is known to satisfy (p^ - vr^) > {K*^ - K^) ~ {D*^ - D^) ~ {B*^ - B^). The 
available lattice data fail to reproduce the approximate plateau exhibited at 
large quark masses. Since we have to reproduce a jlijl25{ri2) interaction, the 
hyperfine splitting is a natural candidate for the improvement. 

To show the trend in the splitting, we compute the average derivative 
< d{m? — m1)/dm1 > for the various cases of interest. This is done via the 
linear fit 

(ml - ml) = Sml + Cml{0) (10) 

(note that both S and C are dimensionless) The summary of these results is 
given in Table §. C should be consistent with 1 when both ml and m^ are 
linear in mg. In this respect, it is interesting to note that C is roughly con- 
sistent with 1 also when we the quark masses are large. The slope decreases 
with /3 both for the standard and for the improved Wilson actions. Again, 
the results at P = 5.7, improved, are consistent with the ones at /5 = 6.3 
, standard Wilson, while the improved data ai (3 = 6.0 have definitively 
the smaller slope, closer to the experimental results (however, still inconsis- 



11 





6.0 5.7 


6.3 W 6.0 W 5.7 W 


s 

K 


-0.027(9) -0.088(23) 
0.92(6) 0.88(5) 


-0.078(11) -0.173(18) -0.207(17) 
0.91(2) 1.01(2) 0.97(2) 



Table 9: Results of the fits for (m^ - ml) = Sml + Cml{0) 

tent with them). Anyway, since the slope is a (slowly varying) function of the 
quark mass, a more detailed comparison is done by simply superimposing the 
results in units of the extrapolated p mass. This is done in Fig. 5. We show 
in Fig. 5a the improved data alone, and the collection of results in Fig.5b. 

The results for the Wilson action at /? = 5.7 are more steeper than the 
other ones in the overlapping region of vr masses (0.2 <~ m^ < 2.0), while at 
f3 = 6.0 and 6.3 the residual difference in the local slope is very small (from 
the plot it is clear that the difference in the average slopes comes mostly 
from the contribution at large masses at f3 = 6.3, and from the one at small 
masses at /? = 6.0). Turning to the improved action, the results at /5 = 5.7 
closely follow the ones at /3 = 6.3, Wilson. The slope of the data at /3 = 6.0, 
which can be compared to the one of /3 = 6.3, Wilson up to m^ ~ 4., is 
clearly more flat, so closer to the experimental data. 

The standard Wilson results seem to have a very slow evolution with 
p. In addition the results at /? = 6.2 reported in |T^ which explore up to 
(7r/p(0))^ ~ 30 do not show any reversal in the decreasing trend , while the 
UKQCD collaboration reported no improvement in the hyperflne splitting at 
the same f3 value with the Clover action [§ (this in the small mass range). 
On the contrary, the data for the improved action from (3 = 5.7 to /? = 6.0 
are still appreciably moving in the right direction, and it is not unreasonable 
to expect good results at a feasible P value (i.e. 6.3 — 6.5). Interesting results 
about the splittings are the ones obtained by the FNAL group 0, which use 
the action ||^. They observe that the spin splittings are very sensitive to the 
parameter c in the 0{a) correction term, —i/2cil)'L^j_yF^yil) which contributes 
mainly an additional magnetic interaction to the quarks. 

Before turning to the baryons, we discuss the results for the meson decay 



constants which are given in Table ^ As explained in [|T^], we obtain /, 



from the local amplitudes in the vf channel whose computation has been 
described in Sect. 3 above. The results are also plotted in Fig.5a, and again 
we compare with the Wilson ones in Fig. 5b. 
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/3 = 5.7 




k 




0.186 0.190 0.194 0.198 


F^MeV) 

Jtt 

J p 


243.5( 55) 219.3( 54) 191. 5( 56) 160. 1( 69) 
0.712( 16) 0.708( 17) 0.721( 21) 0.759( 32) 

0.281( 11) 0.309( 13) 0.341( 19) 0.383( 47) 


/5 = 6.0 




k 




0.176 0.180 0.184 0.188 


F^{MeV) 
U 

Jtt 

J p 


325.9( 35) 293.6( 38) 256.3( 39) 206.4( 39) 
0.669( 7) 0.614( 7) 0.561( 8) 0.521( 9) 

0.182( 4) 0.205( 6) 0.235( 9) 0.276( 20) 



Table 10: Results for the pion and p decay constants. 



The results obtained by making use of the Wilson action aX (3 = 6.0 and 
(3 = 6.3 are in moderate disagreement, while we observe a perfect coincidence 
between the results of the improved action ai (3 = 5.7 and (3 = 6.0 with the 
ones at /? = 6.3, Wilson. This suggests that the data at /5 = 6.3 Wilson or, 
equivalently, 5.7, 6.0 improved, are the asymptotic quenched ones for /^. 

Note also that we convert to MeV by using the extrapolated value for 
the p mass. The /^ data nicely extrapolate at 132 MeV thus demonstrating 
the consistence between the two scales induced by the p mass and /^. (We 
have to say that since we rely on the perturbative evaluation of Z^ 0, it 
is possible, even if unlikely, that this good result is a coincidence. It is clear 
that a non-perturbative computation of the multiplicative renormalization 
constant would be welcome. ) 

We can get an estimate for /^r also from the data in the vr channel, modulo 
a constant. We denote this estimate (/^), and we quote its ratio with the 
true /^, which is remarkably stable. 

fp is more delicate since with Wilson action there is a strong discrepancy 



Even if in 



between the perturbative and lattice evaluation for Zy 
this case the situation is supposed to be better, the perturbative results for 
Zy |2^, hence for the lattice estimate of fp, should be considered with some 
extra care. 

The results for baryons are shown in Figs. 6 in the form of Ape invariant 
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mass plot. The plots we believe are self-explanatory, and demonstrate the 
consistency of the results for the improved action at /3 = 5.7 and j3 = 6.0 
with the ones for the Wilson action at f3 = 6.3 and /3 = 6.0, while the Wilson 
results at /3 = 5.7 exhibit clear scaling violations. 

A final comment concerns the Proton-A splitting. In this case we can 
compare data from different lattices by normalizing with the running value 
of the p mass, like in the Ape plot , thus avoiding the systematics connected 
with the chiral extrapolation. Looking at Fig. 7, we see that the A and the 
proton masses are clearly well resolved, and the results, in the largish errors, 
agree with the ones obtained on the other lattices. ( As discussed in ref. [§] 
the data at /5 = 6.3 are possibly biased from the opposite parity partner.) 
All the other data are interpolating almost linearly between the and the 
infinite mass limit. 
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Figure Captions. 

1 Effective proton mass estimator as a function of tlie time distance at /9 = 

5.7 for tlie different quark masses. 24^ x 32, squares. 9'^ x 18, circles. 

2 a) Pion Green functions for tlie four different quark masses aX f3 = 6.0. 

The results of two (one) particle fits are shown as dashed (dotted) lines. 
b) as a), for the proton. 

3 TT (squares, dashed lines) and vr (diamonds, dotted lines) chiral extrapola- 

tions at /9 = 5.7 (left) and P = 6.0 (right). The lines are the results of 
the fits. 

4 a) The splitting p^ — tt^ plotted against vr^. b) as a), with included the 

Wilson results. Everything in unit of the extrapolated p mass. 

5 a) /tt in MeV versus vr^ in unit of the extrapolated p mass, b) as a), with 

included the Wilson results. 

6 a) Ape plot : nip/ nip vs iri^/ni^ b) as a), with Wilson results. 

7 (rriA — nip) /nip vs iri^/ni^ for the improved and standard Wilson action. 
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